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Abstract 

We consider in this paper a plasma subject to a strong deterministic magnetic field 
and we investigate the effect on this plasma of a stochastic electric field. We show that 
the limit behavior, which corresponds to the transfer of energy from the electric wave 
to the particles (Landau phenomena), is described by a Spherical Harmonics Expansion 
(SHE) model. 



1 Introduction 

This paper is concerned with the effect of a stochastic electric field on a plasma subject to a 
strong magnetic field. This is motivated by the study of the electric turbulence in a fusion 
machine as a Tokamak. Tokamaks are used to confine high energy plasmas in order to obtain 
the conditions needed for nuclear fusion reactions to take place. The plasma evolves in a 
toroidal reactor and is confined in the heart of the torus by the the mean of a strong magnetic 
field. A classical approximation is to suppose the ions to be at rest. Then only the electrons 
are moving. Another classical approximation argument in this type of study is the following: 
we are here interested only in interactions of particles over short distances of the order of the 
Larmor radius, moreover we suppose that at this scale the curvature of the magnetic field- 
lines can be neglected and that the plasma can be considered to be homogeneous along these 
field-lines. Thus we can restrict ourselves to a bidimensional problem. In this approach, the 
Vlasov equation describing the evolution of the repartition function / of the electrons is: 

+ v ■ V x f) + q (Bv 1 - + W turb (i, *)) • V„/ = 0, (1) 

m stands for the electron's mass, q its electric charge, / the distribution function on 
(t,x,v) E K + x M 2 x R 2 , with t the time variable, x the space variable and v the ve- 
locity variable. B is the (constant) norm of the transverse magnetic field, W turb is the 
turbulent electric field, v 1 - is the velocity vector after a rotation of ir/2. We denote 

m 

the cyclotronic frequency, and we want to study the effect of Vy turb in the limit e going to 
zero. In the deterministic case, the limits of related problems have been studied by several 
authors. In the case of the Vlasov-Poisson system (when the electric field is coupled with 
the density J fdv), the limit has been studied, even in the 3D framework, by Frenod and 
Sonnendrucker 0. In the 2D framework, using a slow time scale adapted to the problem, 
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the convergence of the averaged motion, J R2 f(t, x, v) dv to the 2D Euler system of equations 
has been performed simultaneously by Brenier [2], Golse and Saint-Raymond 8 and Frenod 
and Sonnendrucker j^j. A general result in 3D taking into account the two effects has been 
performed by Saint-Raymond in |12| . 

In our case we neglect the Poisson non linear effect, concentrating on the stochastic 
behavior of the equation. Hamiltonian chaos method suggest that the modes of the turbulent 
electric field interacting with the electrons are those having a frequency of u n — 2irne with 
n an integer. This is roughly speaking the Landau resonance. Then the quasi-linear theory, 
(see Garbet jjj) predicts a diffusive behavior with respect to the velocity variable. The 
diffusion coefficient obtained by this method being constant, it can not take into account 
the abnormal diffusion phenomena. In this paper we are interested in turbulent electric fields 
whose spectrum is spread around the Landau frequency and whose spatial fluctuations are 
of the scale of the Larmor radius (of order e). We will show that the limit system is then 
governed by the following equation: 

d t p - d e (a(e)d eP ) =0 (2) 

where e = |i>| 2 /2, p is the average of / over a sphere \v\ 2 = 2e and the diffusion parame- 
ter a(e) is an explicit function of the correlation of V tulh , the turbulent electric potential, 
and of the energy, thus allowing abnormal diffusion. This diffusion parameter is undi- 
mensionally defined by (JSJ). This equation is similar to the so-called Spherical Harmonics 
Expansion (SHE) model in high field limit modeling microelectronics semiconductor devices 
(see P.Degond [S] or Ben Abdallah, Degond, Markowich and Schmeiser pQ). It describes the 
Landau phenomena of transfer of energy from the electric wave to the particles. This work 
uses the techniques introduced by Poupaud and Vasseur ^T] to derive diffusive equation 
from transport in random media. This method works directly on the equation and, for this 
reason, is different from the method used in previous works (see Kesten and Papanicolaou 
[H] and Fannjiang, Ryzhik and Papanicolaou [I]). The paper is organized as follows: 
the precise result is stated in Section|3 In Section|3we show how we can compute explicitly 
the diffusion coefficients. Finally we give the proof of the theorem in Section 0] 

2 Results 

In the remainder of the paper we fix n and we denote 

W e (t, x) = ViW turb (27me£, ex) (3) 
for the stochastic potential. Equation Q takes then the following undimensional form: 

| + „. VJ , + (!i + i=vn^,^v,/, = a (4) 
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we denote E the expectation value of any variable and make the following assumptions on 
the electrostatic potential: 

(HI) V £ £ L 00 (R+;^ 3 ' 00 (M 2 )) and N(e) := E ((HV^Hio.^.-)) 3 ) < oo, 

(H2) EV e (t, x) = 0, for all t G M+, i£l 2 , 
(H3) V e (t 7 x), V e (s, y) are uncorrelated as soon as \t — s\ > 1, 
(if 4) E(F e (t, x)V c (s, y)) = A(t -s,x-y)+ g*(t, s, x, y), 
with: 

^-A £ L°°(R x R 2 ), for a £ N 2 , \a\ < 3, 
ax a 

II |V 2 ,, 5 e | + \Vl,y, y g e \ |U» (C R + ) 3x »«) ^ o, 

where V 2 a g e is the matrix (d Xi d yj g e ). . and V^ yj!/ 3 e is {d Xi d yj d yk g £ ).^ fc . 

Those assumptions are the same than in JJ. Hypothesis (HI) is an assumption on 
the regularity of V e for e fixed. Indeed the norm N(e) can go to infinity when e goes to 
0. Hypothesis (H2) fixes the averaged potential at which is not restrictive. In view of 
(J3J, Hypothesis (H3) determines 27me as the decorrelated lapse of time for the turbulent 
electric field. Namely, it is the bigger lapse of time t—s such that the electric fields at time 
t and at time s can be dependent on each other. Finally Hypothesis (H4), which is very 
classical, can be seen as an homogeneity property which takes place at the local scale e, 
since a quadratic quantity which depends on four variables (t,x,s,y), at the limit, depends 
only on two variables (t — s,x — y). 

We denote R s (v) the rotation of angle s with center of v. We consider the angular 
average of A: 

1 f 27T 

A(t, x) = — / A(t, Rgx) d0. 
27r Jo 

We then have the following result: 

Theorem 2.1 Let V e be a stochastic potential satisfying assumptions (H) and independent 
of the initial data /f G L 2 (R 4 ). Let a(e) be the function defined by: 

a(e) = — 2 J o (-a 2 i)(--^,2^|sin£|) d , (5) 

This function is non negative. Assume that there is a constant Cq such that H/f ||l 2 (k 4 ) < Co 
and: 

e(l + N(ef)^0. (6) 
Let p e be the gyro-average of f e defined by : 

1 [ 27r 

Pe(t,x,e) = — f e (t,x,Rgv)de, (7) 
27T Jo 

for every v such that H 2 /2 = e. Then up to extraction of a subsequence, E/ e ° converges 
weakly in L 2 (R 4 ) to a function f £ L 2 (R 4 ), E/ e converges weakly in L 2 to a function 
f £ L°°(R+,L 2 (R 4 )), Ep e converges in C°([0, T], L 2 (R 2 )-w) for all T > toward a function 
p € L°°(M+; i 2 (R 2 x R+)) n C*° (R+ ; L 2 (R 2 x R+) - W) with p(t = 0, x, e) = L, 2=2e f° dv. 
This function p is solution to: 

d t p~ d e (a(e)d e p) = t > 0, x £ R 2 , e G R+, (8) 

in the distribution sense. Finally E/ e converges weakly in L 2 (R + x R 4 ) toward p(t, x, |u| 2 /2). 
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Remark: depending on the regularity of the function a(e) the solution of the Cauchy problem 
for equation JSJ may be unique. In this case, the whole sequence p e converges to p the unique 
solution of 1(5} . 



3 Explicit computation of the diffusion parameter 

In order to explicit the behavior of a(e) we must need the correlation function A. We assume 
that it follows a "Richardson-like" law 

A(t,x) = f(t)\x\ a . 

Then we have 

a(e) = -— 2 j Q 8lA(^2V2e\ sin -|) ds 

2 fa e a/2 ,+°c _ 



■J n^)i-o 5 r* 



= Ke°' 2 . 

A necessary condition for a function of the form 

p(t,e)=7(i)po(e/f") 
to be an auto-similar solution of (0 is that 

we then have an abnormal diffusion in 

2 

e = t 4 --" . 

For example if a = 4/3 we find a(e) = Ke 2 ^ 3 and an abnormal diffusion in e = i 3 / 4 . 

4 Proof of the result 

We denote 5(K 4 ) the Schwartz space and <S'(IR 4 ) its dual. We denote (•; •) the duality 
brackets between those two spaces. We recall that L 2 (R 4 ) C 5'(K. 4 ) and by extension we 
will denote (•; •) as well for the scalar product on i 2 (E 4 ). For every linear operator P on 
<S(R 4 ) we will denote in the same way P its extension on <S'(R 4 ) defined for every ip <G 5(M 4 ) 
by: 

{Pi/>;n) = (^r,P*v), ??e 5(R 4 ). 

Finally we will say that ip n € 5'(M 4 ) converges to ip E S'(R 4 ) in S'(R 4 ) if for every rj € <S(R 4 ), 
(i/j n ; rj) converges to (tp; if). (This is the weak convergence for <S'(R 4 ).) 
Let us rewrite equation in the following way: 

Btfe + Cf e + ^ = -B\U 

/e|t=0 = fe 
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where C,B,8^ are linear operators on <S(R 4 ) defined by: 



C = v ■ V x 
B = v 1 - ■ V v 

y/e 2itne e 

Notice that C and B are deterministic and non dependent on e nor on t unlike 0\. 

We introduce the projection operator J defined on <S(R 4 ) which averages the values of 
the function on the spheres M 2 /2 = e. Namely, for r\ £ 6>(M 4 ): 

1 i-27T 



Jrj(x,v) = — / r)(x,R e v)d9. 
27T Jo 

We call J the "gyroaverage operator" . This operator is self adjoint for the L 2 scalar product. 
Applying the projection operator J on O and taking its expectation value leads to: 

a t E( J/ e ) + E(JC/ e ) + E(JjB/6) = -E( J0 f e / £ ). 

e 

We first study some properties of the operators in order to pass to the limit in the left hand 
side of this equation. Then we investigate the limit of E(J6|/ e ) following the procedure of 
[TT] . Finally we derive the SHE equation giving the explicit form of the diffusion coefficient 
a(e). 



4.1 Properties of the operators 

We have the following properties on the operators C, B, d\ and J: 

Lemma 4.1 Operators C ,B and Q\ are skew adjoint for the L 2 scalar product. Operators 
C , B , and J commute with the expectation operator E. The operator J is the restriction on 
<S(R 4 ) of the orthogonal projector o/L 2 (M 4 ) into KciB. In particular: 

ll^ll£(L 2 (R 4 )) = 1) 
J 2 = J, 

KeiB = Im J. 



In addition: 

JB = JCJ = 0. 

Proof. 

-Operators C,B and Q\ can be rewritten as b.D, where D is a gradient operator and b a 
regular function verifying D ■ b — 0. For every functions 771 , T72 G iS(R. 4 ) we have: 

(b-D Vl : V2 ) = -( m ;D(br, 2 )) 
= — (771 ; 6 - Dr) 2 ). 

Hence they are skew adjoint operators for the L 2 scalar product. 

-The operator J is clearly the L 2 projection on L 2 functions which depends only on |w| 2 /2 
with respect to v. In polar coordinates v = (rs'm0,r cos8), we have B = d/dO. So J is the 
projection on Ker_B. 
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-Operator J is the projector on Ker_B, hence BJ — 0. Since B is skew adjoint and J is self 
adjoint, we have (BJ)* = —JB = 0. 

-Let us fix 77 6 <S(R 4 ). We denote: Jr)(x, v) = p v (x, ^J-). Hence 

/ Id 2 1 f 27r \ 

jcj v (x,v) = v x • U,o, 1 ~y)^ J o RevdOj = 0. 

Finally JCJ = 0. 

-Since C, B and J are linear and deterministic, they commute with E. □. 

From those properties we deduce the following proposition: 
Proposition 4.2 For every e and every t £ R + we have: 

||/e(*)IU 2 (M 4 ) = ll/elU 2 (M*)- 

There exists a function f° 6 i 2 (M 4 ) and a function f £ L°°(R + ; L 2 (M 4 )) such that, up to a 
subsequence, E/ e ° converges weakly in L 2 (M 4 ) to f°, E/ e converges weakly in L 2 ([0, T] x R 4 ) 
to / for every T > 0. -For every t > 0, / verifies J fit) = f(t). The function JE/ e is 
solution to: 

d t JEf e + E(J9tf e )=w e 

(10) 

JE/ e | t=0 = JE/ e °, 
where w e converges to in S' . 

Proof. Since C, B and 0f are skew adjoint operators, we have: 

d t (f £ (t);Ut))=0, 

which gives the first equality. By weak compactness there exists two functions f° £ i 2 (R 4 ) 
and / £ L°°(E + ; L 2 (R 4 )) such that, up to a subsequence, E/ e ° converges weakly in L 2 (R 4 ) to 
/ , E/ e converges weakly in L 2 ([0,T] x R 4 ) to / for every T > 0. Notice that tQ\ converges 
to in S'(R 4 ). Multiplying Equation © by e, taking its expectation value, and letting e go 
to 0, we find: 

Bf(t) = on R+, 

since B and E commute. Thanks to Lemma 14.11 f(t) £ ImJ, and since J 2 = J, we have 
Jf(t) = f(t) f° r almost every t > 0. Since JB — 0, applying the operator J on equation © 
and taking its expectation value gives: 

d t JEf e + -E(je e J e )=w e , 

with w e = — EJC/ e . This converges in S 1 to —JCJ = —JCJf = 0, thanks to Lemma f4. II 
□ 

Hence we are now concerned by the limit in S' of E(J0f/ 6 ). 



4.2 Computation of E( J0 e J e ) 

Let us denote SI t £ R the group on (S(R 4 ) generated by the operator C + B/e. Namely, 
for every h £ iS, S^h is the unique solution on R to: 

d tg + Cg+^ = (u) 
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The operator S% can be explicitly given by: 

Sth{x,v) = h{T e {t){x 1 v)) 1 

where 

T e (t)(x, v) = (x + ev 1 - - eR^ t/e v x ,R_ t / e v). 

The function T e (t)(x, v) gives the position at — i of the particle being in x with speed v at 
time and moving at constant speed |i>| on a circle of radius e\v\. In particular S% is 2ire 
periodic. Notice that the adjoint of S% is St t - 

Following the procedure of we use a 2 times iterated Duhamel formula. The first 
iteration gives: 

/>27rne 

f e (t) = S^ n J e {t - 27me) - / {Siei^S^SlfS - a) da 

Jo 

and then we write the Duhamel formula for the f e (t — a) in the integral, and this yields: 

/>27rne 

/ (S^U a St a )S^ n Mt-47rne)da 
Jo 

/ / (Sl6l_^_ a ){Sl +a ei_ a _ s S t _ s _ a )St +a Ut-a- S )dsd<j. 
Jo Jo 

We obtain: 

E(J6$f e (t)) = JE(^5| 7rne / e (t-2 7 rne)) 

/•27rne 

-/ E(Jdl(S^t_ a SL a )S^ n J e (t-4wne)) dn + r\ (12) 
Jo 



f 2-nne pA-nne — a 



with 



r; 



/ / JE {eiiSiet^St^SUJU^SU-^Sl+JS -a-s)) ds da. 

Jo Jo 



The function / e ° is independent of the operators d\, ( £ M. In particular, in view of the 
assumption (H3), d\ and f e (t — s) are independent as soon as t > s + 2ime and s > 0. 

Combining this fact with (H2), Equation l|12|l becomes for t > Airne 

- / E(J^(^USlJ)E(^„ e / £ (i - Airne)) da + r t £ , 
Jo 

p2-nne 

E(J0 t e / e ) = - E{J6l{Siei_ a St a ))EfS)da + rl + el, (13) 

Jo 

with e\ = - E(jet(S e Jl_ a St a ))(ESt„ n J e (t-4:7rne)-Ef e (t))da. 
Jo 

Since S% is 27re periodic, S'| wne / e (i — Anne) — f e (t — Anne). We have: 
Siet_ s S*_ s = ±=(SIE%t- s ))-Dl 
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where we denote 

=W e (-^-,-) ) 
2ime e 

and we define the differential operator D e s by 

Dl = R- s/e V v + eR- s / e V^ - eV^. 

Note that Dl is skew adjoint. Let us introduce the operator L\ on <S(R 4 ) (extended on 
S'(R 4 )) defined for every 77 £ S(R 4 ) by: 



li v = - E(et(s^i_ a si a )) v da. 

Jo 



We can gather those results in the following way: 
Lemma 4.3 We have the following equality: 

E(J6if c ) = JLtEf e + rl + et, 
where the operator L\ is defined for every 7] £ <S(M 4 ) by: 

L\n{x, v) = — V„ ■ (E(S e a E f -(t - a) ® £ e (t)) ■ D^(a;, «)) da, (14) 

e Jo 

and </ie remainders are defined by: 

el = JLl(Ef e (t-*nne)-Ef e (t)), 

rt = —}= / (E<(t) ■ V v (S*E*(t - a) ■ D%{Sl +a E\t - a - a) 

e v e Jo Jo 

■Dl +a (St +a f e (t-a-s))))) dsda. 
We can now show the following lemma: 

Lemma 4.4 For every r\ £ <S(R 4 ), the remainder r\ verifies: 

\(rhv)\ <C( V )V~eN(e), 
and {Ll)*rj converges in L 2 (R 4 ) to: 

J R- S V V ■ ( V L^(^ V - R^v^V^ ds. 



Proof. We have: 



( L t)*V = — Dl.(E(S^(t-a)®E*(t))-V vV )do- 
6 Jo 

t>2irn 

= - D%- (V(Sl a E*(t - ea) ® E*{t)) ■ W v ) V da. 

Jo 

But thanks to the definition to T e (s), E e and Hypothesis (H3), the term 

E(S^ E e (t - ea) ® E e (t)) = E^^,^ + v ± - R-^) ® Vff-^,^)) 

2wne 2ime 
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\Dl s ,Dl a E*{t,x)\ < ||VL^ e |U~(K+xR*) 
\Dl,[E^T(e S ))}\ < Ce\V x E*\(T(es)) < C\\ VLnU-m+xM 2 ))- 



converges strongly to {~W 2 xx A)(-a/(2TTn),v 1 - - R- a v^) in L°°((M+) 2 ; VK 1 ' 00 ^ 2 )). Hence 
thanks to the definition of D\, (Lf)*r) converges strongly to: 

J R- S V V ■ (vl x A(^,v x - RsV^Vv'nj ds 

in L°°(R+ x M 4 ). We recall that we have D\ s = i?_ s V„ + ei?- s V^ - eV^ thus 

ll-C , £s'^ > IU 2 (R+xR 4 ) < CII'i'llwL^R+xR 4 ), 
||-De S '-Des ( I ) ||L 2 (R+xR 4 ) < C||^||w 2 > 2 (K+xK*)) 

hence 

Kr^}! < C(n)V~e\\nw^\\fe\\^ 

snp{E(|^(T( eS '))l (\E e (T(es))\ + \Dl s ,[^(T(es))}\) 

s,s' 

(|£ £ (M)I + \D\ s E\t,x)\ + \D\ s ,E%t,x)\ + \Dl s ,D\ s E*{t,x)\))} . 
We then use the following bounds: (recall that E e (t, x) = VF^^, f ).) 

\E e \ < ||V2;y e || L oo( R + xR 2)) 

\Dl s E\t,x)\ < ||V 2 ^ £ ||l- (R+ xr 2 )) 

=-(R+xR 2 )) 

Hln2 Trell 

= (R+xR 2 )) 

Then Hypothesis (HI) ensures that 

l(^}| ^c^ll^lU 2 ^)!!^!^, 2 , 

which ends the proof of the lemma. □ 
We can now state the following proposition: 

Proposition 4.5 Assume that e(N(e)) 2 converges to when e goes to 0. Then the conver- 
gence (up to a subsequence) of JE/ C to f holds in C°(R + ; L 2 (R 4 ) — w), and f is solution 
to: 

d t f + JL a t Jf = 0, (15) 
where the operator L\ is defined for every r\ £ 5(K 4 ) by: 

(L Q t y V = J™ R_ S V V ■ (v 2 ^(^-l,^ - R_,v x )V v t?) ds. 

Proof. Thanks to the previous lemma, for every test function n E 6>(R 4 ): 

1(^77)1^0, 

and (LI)*!] converges strongly in L 2 (R 4 ) to (L^)*rj. But, thanks to Proposition 14.21 f e 
converges weakly to / in L 2 — w. So 

(JLtEf e] r ] ) = (Ef e ;(Lt)*Jr)) 

converges to: 

(f;(L t yj V ) = {JL° t f;r 1 }. 
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The function f e (t) — SWne/e^ — 2irne) converges to in L 2 — w as well. So e\ converges 
to in S'(R 4 ). Passing to the limit in equation (|10fl gives equation (|15fl . This shows that 
dtfe is uniformly bounded in time in a negative Sobolev space. Hence / e converges to / in 
the space of continuous function in time with values in this Sobolev space. Finally since f e 
is bounded in L°°(R+; L 2 (R 4 )), the convergence holds in C°([0, T]; L 2 (R 4 ) - w) for every 
T>0. □ 

4.3 Convergence to the SHE model 

Since J f = /, we can introduce the gyroaverage function defined by: 

p(t,x,e) = f{t,x,v), 

for every v such that 2e = \v\ 2 . This subsection is devoted to the proof of the following 
lemma: 

Lemma 4.6 The function p lies in C°(R+; L 2 (R 2 x M+) - w) (1 i°°(M+; i 2 (M 2 x R+)). It 
is solution to: 

d t p - d e (a(e)d e p) = 
1 f 2n 

Pl*=o = ^ J f(t,x,Rev)d0 
where the diffusion parameter is defined by: 

r2irn f>2iT 

a(e)= / / R e v (-Wl x A)(-— ,R e v ± - R- s+e v x ) ■ R- s R e vdsd6, 

Jo Jo 27rn 
for every v such that e = \v\ 2 /2. 

Proof. Let us first compute the operator JL®J. Let 771 , 771 be two test functions in 6>(R 4 ). 
We have: 

(m;JL t J V 2) = ((L° t yj m ;j m ) 

= / / V^Jr?i(-V^A)(- — ,v ± - R- s v ± )R_ s W v JT] 2 dsdxdv. 
Jr* Jo 2lTn 

Let us denote p Vi for i = 1, 2 the functions defined by: 

Id 2 

Using polar coordinates and noticing that dv = d6 de we find: 

p poo rlun /■27T 

(m;J L t J m) = d e p ni {x,e)d e p ri2 (x 1 e) I I R e e 

J IR^ « 

72 .w-a 



Jo 



R 2 JO 



(— V xa .A)(- — , -R- s+e+ iLe + Re+^e) ■ R^ s+g eds d6 dedx 
Zirn 2 2 

d e p Vl (x, e)d e p m (x, e)a(e) de dx, 



where e = (\/2e, 0). Hence for every test function p^, let us multiply it by Equation (|15Jl 
and integrate with respect to x, v. Since de dO = dv we find: 

/> />00 f- />OG 

&t I / p(t, x, e)p ri (x > e) dxde = / / p(t, x, e)d e (a(e)d e p ri (x, e)) dedx. 
Jr 2 Jo Jr 2 Jo 
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This, with Proposition ^. 51 gives the desired result. □ 
Remark: we have a family of equations parametrized by x G M 2 , and the solutions of two 
equations at two distinct x do not interact. 



4.4 Explicit computation of the diffusion coefficient 

We derive in the following a suitable form to the diffusion coefficient a(e). We will show, 
in particular, that a(e) is non negative. From (H4) the correlation function A(t, x) is even 
with respect to t and x. This with {HZ) gives: 

Lemma 4.7 The correlation function A satisfies: 

SuppA c [-27rn,27m] x M 2 , 
V x A(0,0) =0, 
S A(O,O) = 0. 

This last subsection is devoted to the following proposition. Theorem 12 . II follows from this 
proposition, Proposition 14 . 21 and Proposition 



Proposition 4.8 Let us denote 



I /-27T 



A(t,x) = — / A(R e x,t)d6. 

ZTT 



Then <z(e) is non negative and equal to 

i r , ^ , -s 



/ {-d 2 tt A)(— ,2V^Vl-coss) ds. 
Jo 27m 



27rn 2 

Proof. Thanks to Lemma f4. 61 and lemma l4~7l we have 

poo />2tt 



Since 



f f s 
o(e) = / / R g v ■ (-V 2 A)(-- — , Rev 1 - - R- s R 6 v^) ■ R- s R e v ds d9. 

Js=0 JO 2nn 



-VL^(-^-,^ - R-sV X ) ■ R- S v 
Zirn 

-^V^^--^-,^ - R. S V^) + d 3 (V x A(--^, V 1 - - R-sV^ j), 

Zirn Zirn Zirn 



— / RgvV x d s A{- — ,R e v ± -Rg- s v ± )dsd8 
znn J s=0 J 2im 

2tt 



we find 

a(e) 

R 8 vW x A(0,0) d9 

roc c'2tt 

27171 

where e = (V2e, 0). Let us do the change of variables s f — 6 — s to get 

a ^^2~nn J J 1 { s < e } Re ^' ^^s A {^^i R e+-K/2e~ R s+n / 2 e) d9 ds. 



1 f°° f 2n s 

■ — / / Rge ■ y x d s A{-- — ,R e+n/2 e- R_ s+g+7I / 2 e) ds d6 
'Ttn J s= o Jo zirn 
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Next we have 

s — 9 1 s — 9 

Ree-V x d s A{-^-—,R e+ T T / 2 e-R s+ T T /2e) = --^^d ss A(-^—, R e+7T / 2 e- R s+v / 2 e) 



-dg ^d s A(^-^-,R g+7t/2 e- i? s+7r / 2 e)| 



Integrating by parts the second term of the RHS gives 

,.2n 



[ [ l{s<6}de \d s A{?—^-,R g+7r/2 e- R s+7r/2 e)\ ds 
Jo Jw I 2 7m J 

/ d 8 l {s < g} d s A(- ,R e+n/2 e- R s+7r/2 e)dsdO 

Jo Js. ^ 7Tn 

p2ir p s 9 

/ / 5 s =gd s A(- ,R e+w / 2 e- R s+n/2 e)dsd9 

Jo Js, 27rn 

r 2jr s — 2-k 

/ 9sA(— ,R v / 2 e - i? s+7r / 2 e) ds 

■/-oo 2 7rn 

/ 9 S A(- — ,R v / 2 e-R s+ ^/ 2 e)ds 
d s A(0,0)ds 



d6 



1 

27m . 
1 

27m 

j_ r 27 r - . s 

27m 

i r' 2Tl - -in 

27m 
1 

27m 
1 

27m 

The first two lines cancel by doing the change of variables s' = s — 2n and the third line 
vanishes thanks to Lemma f4. 71 thus 

a ( e ) = To — ^2 / / 1 {^<o}(-9ss A ){^ > :Re+ir/2e~ R s+7r / 2 e)d6 ds 

(Zirn) j j r zirn 

Doing the change of variables s' = 9 — s gives 

2ir poo 



«(e) = 7^2/ / (-d 2 ss A)(--—,R e+7T/2 {{I-R^ s )e))d9ds 
yzixn) j j Z7m 



w| (-diA)(- — A{ I-R_M) dS - 



Finally 



\(I-R s )e\ = V(|l-coss| 2 +sin 2 s)V2^ 



= v/2(l - coss)V2e 
= 2\fe\J\ — cos s 
= 2-\/2e|sin(s/2)| 

which ends the proof of the second assertion. □ 
Computation of the sign of the diffusion coefficient. 

Here we check the non-negativity of the diffusion coefficient by expressing it in another form. 
Thanks to lemma EfTS! and to hypothesis (H3), (HA) we have 

a ( e ) = o^7/ / Rev{-Vl x A){-—,R B v x -Re- s v^)-R e - s vdsdd 
2N J s=0 J_ 2jrN 27m 
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for all N. Then doing the change of variable s := 9 — s wc find: 

fl ( e ) = ™F / / 1 {e>s}Rev ■ (-V 2 xx A)(-^,R v ± - R s v^) ■ R s vdsd0 
But wc remind that thanks to hypothesis (HA) 



Thus 



a(e) = lim —77 



2ttN 

l {e > s} E ( ^ x V\--—-R s v^) ■ R sV ][y x V%~^—-R e v ± ) ■ R e v] ) dsd9 



2ttN 



2irn 2irn 



Interverting s and 9 we see that we can replace l{e- s >o} by l{ s -e>o} an d finally by adding 
both we obtain: 



2N 

a(e) = lim— E 

V ; e,N AN 



/ V x V e (-—,-R,v L )-R a vda 



which is a positive quantity. □ 
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